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Abstract 

A birth-death lattice gas model about the influence of an environment on the fit- 
ness and concentration evolution of economic entities is analytically examined. The 
model can be mapped onto a high order logistic map. The control parameter is a 
(scalar) "business plan". Conditions are searched for growth and decay processes, 
stable states, upper and lower bounds, bifurcations, periodic and chaotic solutions. 
The evolution equation of the economic population for the best fitted companies 
indicates "microscopic conditions" for cycling. The evolution of a dynamic exponent 
is shown as a function of the business plan parameters. 
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1 Introduction 



Economic cycles (EC) (see http : / / cepa.newschool .edu / het / schools /business .htm) 
have been noticed ca. 1930 [1,2]. They have received some political economy 
foundation [3-8], and are often said to be due to time delayed competition 
between "macroscopic variables", investment, production, profit, credit con- 
ditions, ... basically controlled by the interest rate [7-9]. EC have not always 
been well described or defined [10-12]. Cycles are sometimes confused with 
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oscillations, or even noise, in economic data. The measurement of the "pe- 
riod" of a cycle is even far from what would be expected in standard physical 
measurements. 

Economists admit that EC periods are ill defined [10-12] and do depend on 

what goods, assets, ... are considered. Following averaging procedures one can 
consider some constancy (3-4 years) in some EC like from the 1994 NBER 
table. Sometimes not constant periods are said to be the rule, and their values 
checked for according to various data analysis techniques [13,14]. 

EC have recently received some renewed attention in the econophysics com- 
munity [15-18]. One approach that can be thought of is through the stochastic 
resonance phenomenology in which noise and (or) shocks might lead and/or 
induced a stability state, or some dissipative structure [19], smoothly evolving 
[20], - an oscillating trend. A modern econophysics approach would also bear 
upon self-organization [21,22], a notion already thought of by economists [23]. 

Many analytical and simulation approaches are based on the Kaldor - Kalecki 
model [1,2] or generalizations like in [5,24]. With as many parameters one can 
write coupled Langevin equations for fluxes or rates [25]. The main difficulties 
rests in measuring the coupling parameters. 

A simulation with cellular automata or neural networks might lead to a de- 
scription of EC, but there are many caveat in these approaches since they 
also contain parameters or black boxes which are either ad hoc ingredients or 
simply far from so called " flrst principles" . 

Another modern approach implies a microscopic- like description for such a 
macroeconomic problem. This can be at departure from usual continuous time 
evolution equations but the less so if some algorithm is analytically rewritten 
as here below. 

Our analysis stems from a numerically analyzed model [26] considering the 
evolution of a set of entities (agents or companies) under changing economic, 
or more generally environmental, conditions. A question raised was in fact 
whether the system can have bifurcation points. The model is a lattice gas 
in which mobile entities are described by a scalar variable degree of freedom, 
itself compared to what is called a field [26], - controlling the most proba- 
ble state of the system. The degree of freedom evolution depends also on the 
environment of the entity. The economic evolution of the entities is imposed 
to have a Lamarckian feature as often admitted in economy circles [27]. A 
Verhulst attrition-like term is inserted through the company business plan 
which serves for outlining the birth-death process. A statistical physics ap- 
proach based on the logistic map seems thus suitable for describing so called 
economic-like cycles. A posteriori one might wonder why it has not been used, 
since the logistic map ingredients, in the original Verhulst work, stems indeed 
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from socioeconomic considerations [28] . Recall that a pedagogical example of 
a logistic map application is the interest rate effect on savings account [28] . It 
will be found that such "microscopic" ingredients lead to stable states, bifur- 
cations, periodic and chaotic (turbulent-like) regimes. Therefore it is pointed 
out that so called EC belong to a class of self-organized systems and may have 
simply controllable features. Quite importantly it seems that the macroscopic 
ingredients of usual EC theories can be replaced by more microscopic inputs. 
We are of course aware of the present simplicity of the model. EC are not due 
to a birth and death process only. 



2 ACP model 

Let us precise the basic ingredients of the model [26], as used here. A set 
of entities (called "companies") is originally placed at random on a (square 
symmetry) lattice. Each company i is characterized by a real, scalar degree of 
freedom /j G [0, 1]. This parameter denotes how well a company is adapted 
to the economic environment, itself characterized by a (real) field F G [0, 1]. 
The system is e.g. described by the concentration of companies, q G [0, 1]. The 
concentration is defined in a usual way q = ^^rr^, where NtotU) is the number 
of companies on the lattice at time t and N^ax is the volume (the number of 
sites on the lattice) of the system. The time t is discrete and counted in Monte 
Carlo steps (MCS). The companies may diffuse one lattice spacing at a time. 
Although in the original model [26] the field was space and time dependent, it 
will be taken constant in space and time here. Unlike in [26] , where the system 
was described by the position and state of each individual company, within 
this paper, the system is described through a distribution function N{t, f) - 
the number of companies characterized by / at time t. In this case the total 
number of firms existing in the environment at time t is given by: 



In one MCS a number of companies equal to that at the beginning of the step 
is picked at random and possibly displaced by one lattice site. After each move 
the following changes are allowed, if possible: 

(1) The company may disappear due to "difficult economic conditions" as 
measured by F; the probability of surviving is 



where sel is a positive parameter describing the "selection pressure". 
This means that the largest survival probability is when the enterprise i 




(1) 




(2) 
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satisfies the condition: /j = F; it is a best adapted company; 

(2) companies may merge. If company i has a nearest neighbor j, either with 
a probabihty b both form a new company k on site i while j disappears 
with 

fk = + sign[0.5 - r]^^^, (3) 

where r is a random number, r G [0, 1]; 

(3) or with a probabihty 1 — b, companies survive and create one (or two) 
new firm(s), depending on the available space in the Moore neighborhood 
of the i-company, each with a specific random /. 



3 Mean field approximation 

Averaging over a MCS such that 

< + sign[0.5 - r]i^^ >=< A±A >^ /^(^ + i) 



taking into account the modification in the state distribution caused by (ii), 
i.e. N{t + 1,/) N{t, (/ + F)/2) and considering all possible events in a 
mean field sense the evolution equation of the system can be written as the 
sum of two terms, 

N{t + 1, /) = H^{ct)p{f)N{t, /) + H2{ct)N{t, g{f)), (4) 

where (dropping the t index in Cj) : 

H^(c) = 1 - 6(1 - c)[4c^ + 6c^(l - c) + 4c(l - cf] 

+4(1 -6) (1 -c)c^ (5) 

+2(1 - 6)(1 - c){l - - 4(1 - cfc 
-6(1 - c)^c^ - 4(1 - c)V - (1 - cfc"]}, 
H2{c) = b{l - c)[4c^ + 6c\l - c) + 4c(l - cf] (6) 

and 

f + F 

9if) = (7) 
Eq.(4) can be iterated back until the initial point 

N{t, /) = E a,,,_,p'-^{f)N{0, g%f)), (8) 

1=0 
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1.2 




b 

Fig. 1. Iterations of the concentration as a function of the parameter b starting from 
a low initial concentration 

where each coefficient ai^t-i is the sum of time oriented products of H\{ct) and 
-f^2(ct), i.e. 

O'ht-i = X)-f^a(t)(Cto)-^<T(t-l)(Cti) • • • H„(t-i){Cti)- (9) 
c 

where > t\ > t2 > ■ ■ ■ > ti and a denotes all possible permutations of the 
set {1,1,. ..,1,2, 2,. ..,2}. 

^ V ' ^ V ' 

itimes t-i times 



4 Logistic map 

The evolution equation for the total number of the best adapted companies 
has the form of a logistic map [28,29], with a control parameter h. Therefore 
the model may reveal some periodic or chaotic behavior. The best adapted 
company concentration evolution, starting from a small initial value, is shown 
at fixed 6's in Fig. 1 for different h values. The corresponding Lyapunov ex- 
ponent A is shown in (Fig. 2). Five main regions can be distinguished in Fig. 
2. Only one stable solution (with q > 1) exists for h > 0.45, i.e. due to the 
frequent merging process, there is no overpopulation of the system. The sta- 
ble state is achieved as a dynamic equilibrium between merging and spin-off 
creation process. For other b values some "convergence" is found toward one 
or several states characterized by bifurcations and limit cycles with various 
periodicity down to be — 0.1556, where a chaotic region appears containing 
narrow stability windows. 

In order to better understand the dynamics we calculated 

(1) the time needed before a stable orbit (state) is reached. Fig. 3, and 

(2) a variant of the Lyapunov exponent in order to emphasize both the 
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Fig. 2. Lyapunov (A) (solid line) and the dynamic exponent ( (dashed line) as a 
function of b 

convergence (or divergence) of the trajectory, and a possible oscillatory 
behavior, through the so called real (A') and imaginary (A") part of a 
"generalized" A = ln[{d / dx) f (x)] in conventional notations, but without 
absolute values before taking the logarithm. Fig. 4. 

For b > 0.45 the time required for the system to attain the stable state remains 
(approximately) on the same level until b < 0.65. However for 6^1 the time 
required for the system to reach the stable state increases rapidly (Fig. 3). 
This corresponds to the behavior of the Lyapunov exponent (Fig. 2), which 
is A = for 6 = 1. Especially interesting is the range b e (0.38; 0, 45): damp- 
ing properties are superposed to an oscillating behavior. Oscillations are also 
visible in Fig. 4, where the finite imaginary part of A shows that the system 
may cycle. Notice the value of A" equal to tt, tt/2, 37r/4,... in various stability 
regions, but taking bizarre values in chaotic regions. This is another way to 
find characteristics of stability windows in the "turbulent regimes". 

The process of achieving a stable state has been also studied by assuming a 
power law decay for the c-distance evolution in the form \c{t)\ ~ i^. The plot 
of the C exponent as a function of b is shown in Fig. 2. 



5 Conclusions 

The above results, for a simple model, lead toward interesting conclusions con- 
cerning stable or not economic states. We do not - of course claim that EC 

are due to a merging - spin off microscopic process only. What is intended is 
to search whether a discrete microscopic model has some interesting ingredi- 
ents for some macroeconomic evolution. The evolution of the model depends 
on the three microscopic-like parameters the environmental condition F, the 
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Fig. 3. Time for achieving a stable state as a function of b 
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Fig. 4. Real and imaginary part of a Lyapunov-like exponent 



selection pressure sel and the merging probability b having a physical sense in 
contrast to macroscopic (economic) ones . The influence of the environment F 
and the selection pressure control the time required for the system to achieve 
a stable state. Moreover the environmental condition decides also on the type 
of companies surviving in the system. The first parameter is essentially of 
pohtical and/or global origin. The second can be taken as the inverse of tem- 
perature in thermodynamics, and leads as in other non equilibrium systems 
to consider an intrinsic noise controlling the existence of dissipative structures 
[19,28]. If the selection pressure is "strong" then the "weak" companies are 
removed from the system and the distribution function N{t, f) converges to 
the "best adapted companies" case N{t,F). This is also similar to results 
found in Darwinian-like economic evolutions [30,31]. The last parameter h is 
more subtle, and concerns business plans: it controls the asymptotic concen- 
tration value of companies and the cyclic or chaotic behavior of the system. 
We have constrained 6 to be a constant though a dynamic economic policy 
usually requests a time dependent h. 
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Concentrating on the best adapted companies for small b the system reveals 
chaotic properties, i.e. A > 0. It means that at low merging the system becomes 
unstable. This is an apparently new interesting feature which seems reasonable 
but nevertheless requests economic studies and considerations. The system is 
especially stable at 6 ~ 0.38 where the damping parameter attains its largest 
value. In such a case, the system has only one stable solution, which is achieved 
in the shortest possible time. Therefore a merging philosophy, or politics, 
is a stabilization factor for economic systems. In fine, it should leads to 
a monopolistic situation which is stable and " cycling" as long as the goods of 
such a company are needed, - if there is no spin-off process. 

Finally, the above numbers and findings should truly be taken with caution 
since they pertain to a model studied on a 2D lattice with a given symmetry. 
The period of cycles depends on the economic plan parameter b, but also on the 
Moore neighborhood hereby chosen for the evolution and the information flow 
between companies. A posteriori it is understood why "EC periods", either 
global or for specific activity fields, are poorly defined [3,4,12]. Interestingly 
it appears again that the connectivity of the "lattice" on which an economic 
process takes place is markedly relevant for defining economic oscillations. 
This has already been noticed in the case of financial crashes [22,32]. Whence 
some EC appear to be self-organized systems and have simply controlled (and 
controllable) features, similar to turbulent processes. 
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